Phase Diagram and Fixed-Point Structure of two dimensional J\f = 1 Wess-Zumino 

Models 
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We study the phases and fixed-point structure of two-dimensional supersymmetric Wess-Zumino 
models with one supercharge. Our work is based on the functional renormalization group formulated 
in terms of a manifestly off-shell supersymmetric flow equation for the effective action. Within 
the derivative expansion, we solve the flow of the superpotential also including the anomalous 
dimension of the superfield. The models exhibit a surprisingly rich fixed-point structure with a 
discrete number of fixed-point superpotentials. Each fixed-point superpotential is characterized 
by its number of nodes and by the number of RG relevant directions. In limiting cases, we find 
periodic superpotentials and potentials which confine the fields to a compact target space. The 
maximally IR-attractive fixed point has one relevant direction, the tuning of which distinguishes 
between supersymmetric and broken phases. For the Wess-Zumino model defined near the Gaufiian 
fixed point, we determine the phase diagram and compute the corresponding ground-state masses. 

PACS numbers: 05.10.Cc,12.60.Jv,11.30.Qc 



I. INTRODUCTION 

Supersymmetry has become a well-received guiding 
principle in the construction of particle-physics models 
beyond the standard model. Whereas the resulting phe- 
nomenology of these models is often worked out by per- 
turbative analysis, the required breaking of supersymme- 
try may be of nonperturbative origin; it is therefore typi- 
cally parameterized into the models by a phenomenologi- 
cal reasoning. For a proper understanding of the underly- 
ing dynamical mechanisms of symmetry breaking which 
is often related to collective condensation phenomena, 
powerful and flexible nonperturbative methods specifi- 
cally adapted to supersymmetric theories will eventually 
be needed. 

For many nonperturbative problems in field theory, lat- 
tice formulations and simulations have proven success- 
ful. As supersymmetry intertwines field transformations 
with spacetime translations, discretizing spacetime of- 
ten induces a partial loss of supersymmetry. This prob- 
lem also goes along with the challenge of properly im- 
plementing dynamical fermions on the lattice, currently 
witnessing significant progress 

[IS IIH- As completely 
new territory is entered in these studies, nonperturbative 
continuum methods which can preserve supersymmetry 
manifestly can complement the lattice studies, eventually 
leading to a coherent picture. 

A promising candidate for a nonperturbative method 
is the functional renormalization group (RG) which has 
been successfully applied to a wide range of nonperturba- 
tive problems such as critical phenomena, fermionic sys- 
tems, gauge theories and quantum gravity, see [H, [|| 0, H, 
reviews. A number of conceptual studies of su- 
persymmetric theories has already been performed with 
the functional RG. The delicate point here is, of course, 
the construction and use of a manifestly supersymmetry- 
preserving regulator. For instance, a supersymmetric 
regulator for the four-dimensional Wess-Zumino model 



has been presented in [IJ, [L2| ■ A functional RG formula- 
tion of supersymmetric Yang-Mills theory employing the 
superfield formalism has been given in [13| : for applica- 
tions, see also 0, HEj . Recently, general theories of a 
scalar superfield including the Wess-Zumino model have 
been investigated with a Polchinski-type RG equation in 
[l6l |. yielding a new approach to supersymmetric non- 
renormalization theorems. A Wilsonian effective action 
for the Wess-Zumino model by perturbatively iterating 
the functional RG has been constructed in [17| ■ 

This work is devoted to the two-dimensional Af — 1 
Wess-Zumino model with a general superpotential, ex- 
ploring the model beyond the realm of perturbative ex- 
pansions around zero coupling. This is the simplest quan- 
tum field theoretic and supersymmetric model where 
the nonperturbative dynamical aspects of supersymme- 
try breaking can be studied. The present study details 
and generalizes our results presented in a recent Letter 
[l8| . and builds on our earlier work on supersymmetric 
quantum mechanics, where we have constructed a mani- 
festly supersymmetric functional RG flow for the anhar- 
monic oscillator [l^| ; see also [2(| [2l[ for RG studies of 
supersymmetric quantum mechanics. 

Inspired by Witten's work on the potential breaking of 
supersymmetry in the Wess-Zumino model (22j . pioneer- 
ing nonperturbative lattice studies based on Hamiltonian 
Monte-Carlo methods had early been performed for this 
model by Ranft and Schiller 23] . More recently Beccaria 
and coworkers [U [H| re- investigated the phase diagram 
and the ground-state energy of themodel with similar 
methods. Golterman and Petcher [26[ formulated a lat- 
tice action with a partially realized supersymmetry. An- 
other lattice study of the Wess-Zumino model has been 
performed by Catterall and Karamov [27|. The results 
of these studies give a first glimpse into the properties of 
the phase diagram as will be discussed in Sect. IVIII 

In the present paper, we use the functional RG equa- 
tion for the superpotential W to study the phase struc- 
ture of the Af = 1 Wess-Zumino model in two dimen- 
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sions. In Sect. IIII1 we begin with recalling the off-shell 
and on-shell effective potentials in a one-loop approxima- 
tion and comment on certain flaws of the approximations. 
The main part of this work is concerned with extending 
and applying the manifestly supersymmetric RG tech- 
niques developed in in the context of supersymmet- 
ric quantum mechanics. The manifestly supersymmetric 
flow equation for the effective action is constructed in 
Sect. IIV1 To first order in a derivative expansion of the 
effective action, we solve the RG flow equation for the 
effective superpotential in Sect. [V] 

In the fixed-point analysis performed in Sect. IVII al- 
ready a simple polynomial expansion of the superpoten- 
tial gives access to an infinite number of fixed points with 
an increasing number of RG relevant directions. Beyond 
the polynomial expansion, the flow equation for the full 
superpotential reveals a variety of qualitatively different 
solutions depending on the initial conditions: we find 
periodic, sine-Gordon type solutions as well as sigma- 
model type solutions confining the field values to a fi- 
nite interval. At next-to-leading order, a nonzero anoma- 
lous dimension governs the large-field asymptotics of the 
fixed-point superpotentials, such that a family of regu- 
lar fixed-point solutions arises. This family of superpo- 
tentials shows oscillating behavior for small fields and 
a standard asymptotics for large fields, similar to fixed- 
point potentials for pure bosonic theories in two dimen- 
sions [H, HH ■ As a particularity of these supersymmetric 
models, we identify a new scaling relation between the 
leading critical exponent of the superpotential flow and 
the anomalous dimension. 

Finally, we study the phase diagram of a particular 
Wess-Zumino model defined near the Gaufiian fixed point 
in terms of a quadratic superpotential perturbation in 
Sect. I VIII Depending on the initial values of the control 
parameters of the potential, we observe a quantum phase 
transition from the supersymmetric to the dynami cally 
broken phase. Following the lattice studies 0, [25|, |28| . 
we calculate the critical value of the control parameter 
for the phase transition as a function of the coupling A 
at the cutoff A, but now for all values of the coupling A. 
We also determine the fermionic and bosonic masses in 
both phases. 



The anticommuting 9 parameter in this expansion is a 
constant Majorana spinor. Supersymmetry transforma- 
tions are generated by the supercharges 



a = -4 



(2) 



which anticommute on space-time translations, 
{Q a , Qp} = 20ai3- The transformation rules for 
the component fields are obtained by comparing 
coefficients in (5$ = ie[Q, <&] and read 



5<fi = e7*"0! 6ip = (F + ij*0(f))e, 
54, = e{F - i$<fry* ) , 8F = ie0i/j. 



(3) 



As usual the F term transforms into a total derivative 
such that its space-time integral is invariant under su- 
persymmetry transformations. The supercharges anti- 
commute with the superderivatives 



V = ■£= + iq. 
86 



(4) 



and up to a sign they obey the same anticommutation 
rules as the supercharges, {T> a ,T>p} = —2i$ af} . In ex- 
plicit calculations, one uses Fierz identities which all fol- 
low from 



i>x = -§xV>- sT/iOo/WO - \i*(xi*i>)i 



(5) 



where 7* = «7o7i anti-commutes with the 7^. It is useful 
to keep in mind that for Majorana spinors the fermionic 
bilinears have the symmetry properties 

$X = -X*l>> i>lnX = -Xl^ and ^7*% = Xl*^, (6) 

such that the only Lorcntz- invariant bilinear is rpj*ip 
since iptp = 0. Thus, we choose as Lagrangian density 
£0 f° r the free theory the D term of 

i©$7„Z>* = i^7*V + {h*^)F - ih^d^ ~ 0-yJCo. 

(7) 

In components it has the form 

C = \d^<i> + - jS M ^7"V ~ \F 2 - (8) 



II. WESS-ZUMINO MODEL 

Two-dimensional Wess-Zumino models with one su- 
persymmetry are particular Yukawa models where the 
self-interaction of the scalar field determines the Yukawa 
coupling. In an off-shell formulation they contain a scalar 
field <p, a Majorana spinor field ip and an auxiliary field 
F. To maintain supersymmetry in every step of our cal- 
culations we combine these fields to one real superfield 

$(x, 9) = 0(a) + 6^{x) + h§yJ)F(x). (1) 



In this work, we study a class of interacting theories, 
where the interaction Lagrangian is given by the D term 
of a superpotential Vt / ( < I > ), 

w($) = w(<f>) + e^ipw'((t>) - \h*0Ci. (9) 

In components, C\ has the form 

d = \w"{4>)$i*^ - w'(4>)F. (10) 

The sum of Co and C\ defines the off-shell Lagrangian 
density 



c = ^d^+m 



l_p2 1 1 W " 



2 W"(^^-W'(cb)F, 
(11) 
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which gives rise to an invariant action. As expected for 
a Euclidean model, this action is unbounded from below 
and above. After eliminating the auxiliary field via its 
algebraic equation of motion 

F = -W'(4>), (12) 

we end up with the stable on-shell Lagrangian density 

C = i(^) 2 + tffap + \W 2 ((f>) + \W"{cj>) ^ 7 *V- (13) 

This density is invariant under the nonlinear on-shell su- 
persymmetry transformations 



6(f> = 67* V, H = {ilJ<t> ~ W'{<f>))e, 
= i(ijJ<t)-W' ((/>)). 



(14) 



For a polynomial superpotential W the supersymmet- 
ric Yukawa models defined in Eq. (fTTj) and (fT3")) are per- 
turbatively super-renormalizable. If the leading term in 
the superpotential contains an even power of <p, W = 
c<f> 2n + 0(4> 2n ^ 1 ), supersymmetry cannot be broken by 
quantum corrections. However, if the leading term con- 
tains an odd power, supersymmetry may be broken. 



III. ONE-LOOP PERTURBATION THEORY 

Let us first discuss the effective potential in one-loop 
approximation, which can be set up in both the on-shell 
or the off-shell formulation. It is well-known that the 
one-loop on-shell potential becomes artificially complex 
for nonconvex classical potentials [29T |. This problem is 
avoided in the off-shell one-loop formulation: here, keep- 
ing first the auxiliary field in the one-loop calculation 
and subsequently eliminating it by its quantum equation 
of motion corresponds to a resummation of higher-order 
terms in the on-shell formulation. We expect that this 
potential is a better approximation to the exact effective 
potential as compared to the one-loop on-shell poten- 
tial. Indeed, we find a real and stable effective potential 
based on the off-shell calculation. Similar observations 
can be found in 30] . The corresponding problem in one- 
dimensional supersymmetric systems has been carefully 
analyzed by Bergner (see [3lT | and references therein). 



A. On-shell effective potential 

To calculate the one-loop potential in the on-shell for- 
mulation, we need the fluctuation operators M-p and Mb 
for the fermion and the remaining scalar, respectively. 
For a homogeneous background field <f> playing the role 
of a mean field, these operators read 



M B = p 2 + V" and M F = ic, 



Ml 



(P 2 



+ 1*W" 
W" 2 )l 2 , 



(15) 



where V = is the classical potential for the scalar 

field. In a finite box of size L, the momentum takes 



the values p^ — 2Ttn fJ- /L with integer- valued n^. 
fcrmionic integration yields the Pfaffian of Mp, 



Pf(M F ) = ±y/det(M F ) 



The 



(16) 



We shall assume that the Pfaffian has a fixed sign, such 
that we may replace the Pfaffian by the square root of the 
determinant. In a perturbative approach and in the bro- 
ken phase (with vanishing Witten index), this assump- 
tion is justified. 

For an operator M with eigenvalues p 2 + C 2 the deriva- 
tive of the zeta function Cm(s) = tr(M/ [i 2 )~ s at the ori- 
gin is 



Kt(CLn) 
CLn 



(17) 

with the last sum, which involves the MacDonald func- 
tion, approaching zero exponentially fast with increasing 
box sizes. Thus, in the thermodynamic limit the one-loop 
effective potential in the zeta-function scheme is 

u ® =\ w ' 2 + ^i(Cm p (o)-Ck(o)) 



-w- 



2L 2 
W" 2 

~8^r 



[(1 + Y)ln(l + Y) 
+ X ln(W" / fi) 2 — X] 



(18) 



with X = W'W"'/W" 2 . The energy scale fi is fixed by 
a renormalization condition. If we use a momentum cut- 
off rcgularization instead of the ^-function regularization 
then we obtain the same result with \x replaced by the 
cutoff A. 

The effective potential becomes complex for nonconvex 
classical potentials V = ^W' 2 . To be specific, let us 
choose 



W = A(0 2 - a 2 ) 



V 



A 2 



a 2 ) 2 



(19) 



For negative a 2 , the one- loop effective potential is real. 
For positive a 2 , it becomes complex for small fields 
4> 2 < a 2 /3. For fields slightly bigger than a/v3, the 
potential is real and negative. This signals the failure 
of the approximation since the effective potential must 
be non-negative in a supersymmetric theory. Depend- 
ing on the sign of the renormalized a 2 , we find both a 
supersymmetric phase characterized by a non-vanishing 



expectation value (4>) 



6 min and Uon( 



and 



a phase with ((f>) = and broken supersymmetry, see 
Fig. [T] Here, the renormalization conditions are chosen 
such that the minimum <^ m i n agrees with the renormal- 
ized value a in the supersymmetric phase. In the broken 
phase, we use a simple renormalization condition by fix- 
ing all parameters at the cutoff, fi = A. In both cases, 
the renormalization condition could alternatively be for- 
mulated in terms of Coleman- Weinberg renormalization 
conditions by fixing the curvature of the potential at the 
minimum to the physical mass. 
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B. Off-shell effective potential 

In the off-shell formulation, the fluctuations of both 
scalars (including the auxiliary field) and the fermion are 
taken into account. The off-shell Lagrangian (fTTj) gives 
rise to the fluctuation operators in the background of 
both a d> and an F mean field, 



M H = 



- FW' 
-W" 



-W" 

-1 



and Mf=^ + 7*W' 



(20) 

where the mean fields (as the argument of W" and 
W") and F are assumed to be homogeneous. It follows 
that the ^-function regularized one-loop off-shell poten- 
tial reads 



I W" 2 
C/W = - -F 2 - FW' - + Y) ln(l + Y) 



(21) 



+Yln(W" I [if -Y] 




FIG. 1: The classical potential V and the on-shell and off-shell 



effective potentials Uo^ 
The subscripts + and — 



and Utff in one-loop approximation, 
denote the sign of a 2 . 



IV. SUPERSYMMETRIC RG FLOW 



with Y — —FW"'/W" 2 . With a momentum cutoff regu- 
larization, we obtain the same result with /i denoting the 
cutoff. To eliminate the auxiliary field F, we must solve 
the transcendental gap equation 



-F-W 1 ■ 



W" 
8tt 



In 



(W 



112 



FW h 



(i) 

off 



0, (22) 



For an 



and insert the solution for F back into £/ 
arbitrary </>, the gap equation always has a real solution 
F leading to a real effective potential [32|. Concerning 
supersymmetry breaking, we find the same qualitative 
result as with the on-shell calculation: the sign of the 
renormalized a 2 determines the phase of the system. 

Note that the on-shell and off-shell potentials (fT5|) and 
(|2"Tj) have similar forms. In Fig. [1] the two potentials 
are compared with each other and also with the classical 
potential V. Of course, the same renormalization con- 
ditions are used for the on- and off-shell potentials. For 
positive a 2 , the scale parameter fi has been adjusted such 
that the potentials take their minima at the same value 
0min = a. Note that the off-shell potential contains re- 
summcd contributions from higher order in ft. The reason 
is that the solution F of (f2"2"| contains terms of higher or- 
der in ft and inserting the solution back into the effective 
potential generates terms to all orders of ft in the effec- 
tive potential. By expanding the off-shell action in ft, the 
first-order result agrees again with the complex on-shell 
effective potential. The effective resummation contained 
in the off-shell action is such that the effective potential 
becomes real and non-negative everywhere, in particular, 
at those points where the classical potential is not con- 
vex. For a 2 > 0, both one-loop potentials predict a phase 
with broken TL-i symmetry and unbroken supersymmetry. 
For a 2 < 0, we find a phase with unbroken %2 symmetry 
and broken supersymmetry. 



In this section, we will construct a manifestly super- 
symmetric flow equation in the off-shell formulation. Our 
approach is based on the functional RG formulated in 
terms of a flow equation for the effective average action 
Tk, i.e., the Wetterich equation, [33| 



d k T k = 1 STr 



(23) 



Here, Tk is a scale-dependent effective action; it inter- 
polates between the microscopic or classical action S for 
k — > A, with A being the microscopic UV scale, and the 
full quantum effective action Y = Yk^o-i being the stan- 
dard generating functional for 1PI correlation functions. 
The interpolating scale k denotes an infrared IR regulator 
scale below which all fluctuations with momenta smaller 
than k are suppressed. For k — * 0, all fluctuations are 
taken into account and we arrive at the full solution of 
the quantum theory in terms of the effective action T. 
The Wetterich equation defines an RG trajectory in the 
space of action functionals with the classical action S 
serving as initial condition. 

In Eq. (|23p , we encounter the second functional deriva- 
tive of Tk , 



(if) 



ab 



-A. 



(24) 



where the indices a, b summarize field components, in- 
ternal and Lorentz indices, as well as spacetime or mo- 
mentum coordinates. In the present case, we have 
= (<j),F,ip,ip) where $ is not a superfield, but 
merely a collection of fields. The momentum-dependent 
regulator function Rk in Eq. (|2"3"|) establishes the IR 
suppression of modes below k. In the general case, 
three properties of the regulator Rk(p) are essential: (i) 
Rk(p)\ P 2 /k 2 ^o > which implements the IR regulariza- 
tion, (ii) Rk(p)\k 2 / P 2 ^a — which guarantees that the 
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regulator vanishes for k — > 0, (iii) Rk{p)\k^A^<x> — * oo 
which serves to fix the theory at the classical action in 
the UV. Different functional forms of R& correspond to 
different RG trajectories manifesting the RG scheme de- 
pendence, but the end point I\.^o — ► V remains invari- 
ant. 

The regularization preserves supersymmetry if the reg- 
ulator contribution to the action ASk is supersymmetric, 
see below. As the regulator needs to be quadratic in the 
fields in order to maintain the one-loop structure of the 
flow, a general supersymmetric quadratic form can be 
constructed as a D term of a superfield operator $/v$. 
Here, if is a function of the two invariant and commuting 
operators 2?7*2? and V$D oc (Dj*T>) 2 . Since powers of 
T>^ ie T> boil down to 



2n+l 



2\n-l 



and 



2\n 



(25) 



where d 2 is the standard Laplacian, any invariant and 
quadratic regulator action is the superspace integral of 



1 



*2?(fi(-3 2 ) - 7 *r 2 (-<9 2 ))£>$. 



Expressed in component fields, we find 



AS k = I 



,F)R% 



F 



(26) 



(27) 



the flow equation. In this work, we use the derivative 
expansion which is based on the underlying assumption 
that the fully interacting theory remains sufficiently local 
if formulated in the given set of field variables. In order to 
preserve supersymmetry, we expand the effective action 
in powers of super-covariant derivatives in the off-shell 
formulation. This expansion allows for a systematic and 
unique classification of all possible operators. A trun- 
cation of the effective action to a finite derivative order 
leads to a closed set of equations for the expansion pa- 
rameters. 

In this section, we concentrate on the leading-order 
derivative expansion: the so-called local potential ap- 
proximation. Here, the truncated effective Lagrangian is 
given by Eq. (fTTj) with a scale-dependent superpotential 
Wk, such that the truncated effective action reads 



r k [<t>,F,$,i>]= / tfx^d^t + ^H 



if 2 



+ W£{4>)M-w' k {<l>)F) 



(29) 



The derivation of the flow equation for the superpotential 
parallels the corresponding one for supersymmetric quan- 
tum mechanics given in a previous work [19( . Within the 
approximation of constant mean fields, the second func- 
tional derivative of the effective action plus regulator is 



In momentum space, id^ is replaced by p^ and the oper- 
ators take the explicit form 



Rt 



p 2 r 2 



and 



Rl 



>r 2 + 7*ri, (28) 



where r% = p 2 f\. Comparison with Eq. (|20p reveals that 
r\ plays the role of a momentum-dependent supersym- 
metric mass term, whereas r 2 can be viewed as a de- 
formation of the momentum dependence of the kinetic 
term. 

This choice of the regulator guarantees a supersym- 
metric RG trajectory; i.e., for a supersymmetric initial 
condition I\._>a — > S, the solution to the flow equation 
Tfe will remain manifestly supersymmetric for all k in- 
cluding the endpoint T — T^o- This does not only 
hold for the exact solution, but is also valid for trun- 
cated effective actions, provided the truncation is built 
from supersymmetric field operators. 



LOCAL POTENTIAL APPROXIMATION 



,(2) 



n k — | w/'/„, „/, <> „T 



ei 



(1, of, 



B 



(30) 



where the operators on the diagonal read 



, = (p 2 (l + r 2 ) - FW' k " + ^W^ifa^ -Wfl - n 
-Wg-n -l-r 2 



S = i(l + r 2 )s> + 7*(ri + Wfe)- 



(31) 



The inverse of the operator defined in Eq. (|30|) can be 
written as follows, 



-.(2) 



Rk 







or 


Gl F ) 


A- 1 










K 



A- 1 




x 

B- 1 



(32) 



Various systematic and consistent approximation 
schemes for the construction of Tk can be devised with 



where we have abbreviated 
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K 



A + W^i^B- 1 ^) 



(33) 



r 



In order to verify that Eq. (|32|) is the inverse of Eq. (|30|) . 
the Fierz identity ijj ® ijj = — 57,(^7*?/') is useful. This 
result is inserted into the flow equation (I23|) . which in 
component notation reads 



i Tr (d k R* Gl B ) - i Tr (d fe i? F G FF ) . (34) 



The flow equation for W k (<f>) is obtained by projecting 
both sides of this equation onto the term linear in the 
auxiliary field. This yields 



d 2 p /(l + r 2 )(^' + n) 



W(0)=-WJ" / ^ 



A 2 

P 2 (l + r 2 ) 2 -{W^ + ri ) : 
2A 2 



•' -2), (35) 



where we have introduced A = p 2 (l + r 2 ) 2 + (W k + ri) 2 . 
Integrating with respect to and dropping an irrelevant 
constant leads to 



d k W k {<i>) 



1 



d 2 p (r 2 + l)d k r 1 -(r 1 +W' k '(^))d k r 2 



(2n) 2 



A 



(36) 



This flow equation for the superpotential has exactly the 
same structure as the corresponding flow equation in su- 
persymmetric quantum mechanics [19| . This is not sur- 
prising since supersymmetric quantum mechanics can be 
obtained from this model through dimensional reduction. 

Here we are interested in superpotentials for which the 
map It 3 <f> — > W'(4>) € 1R has winding number zero 
as these potentials allow for dynamical supersymmetry 
breaking. For a polynomial W this is the case if W tends 
asymptotically to an even power, W'{(f>) ~ c<j) 2n . Then, 
the highest power of W"((p) is odd. This implies that the 
mass-like regulator r% does not screen but merely shift 
possible zeroes of WL'. Thus we may set n = without 
spoiling the IR properties of the flow. 

In the present local-potential approximation, the sim- 
ple cutoff function r 2 = (k/\p\ — 1) 9(l—p 2 /k 2 ) turns out 
to be technically very convenient, since the momentum 
integration in Eq. (|36[) can be performed analytically, 



d k W k {<f>) =- 



WW) 



4tt k 2 + Wg(4>) 



(37) 



In order to calculate the bosonic potential V(<j>) = 
^W'((f>) 2 , we only need the derivative of the superpo- 
tential. The corresponding flow is 



d k W' k {4>) = 



W> k " 



(0) 



k k 2 - W' k \<j>) 2 
4tt (k 2 + W^((/>) 2 y 



(38) 



This equation exhibits a particularity for any finite value 
of k, as the sign of the flow depends on whether W k (<fi) 2 
is smaller or larger than k 2 . For large 4>, we generally 
expect both (W k ) 2 and W k " to be large and positive for 
TL 2 symmetric systems. In this case, the flow for large <p 
tends to deplete the height of the potential. Of course, 
for <j) — > 00, we expect the denominator of Eq. (|38p to win 
out over the numerator, such that the flow vanishes at 
large field amplitudes. For small <f>, or, more generally, 
in the vicinity of local or global minima of W' k , there 
can be an inner domain where (W k ) 2 < k 2 . For convex 
potentials W k with W k > 0, the flow is negative here, 
resulting in the tendency to flatten out this inner part of 
the potential W k . As the curvature W k " is related to the 
masses of the excitations, the flow shows a clear tendency 
to small masses if an inner domain with (W k ) 2 < k 2 
exists. As it will turn out later, this is a characteristic 
property of the supersymmetry-broken phase. 

Let us note in passing that the regulator used in the 
present section can lead to artificial divergences at higher 
orders in the derivative expansion, e.g., when a wave 
function renormalization is included. Then a stronger 
regulator in the IR is needed; see, App.lAlfor calculations 
at next-to-leading order in the derivative expansion. 



VI. FIXED-POINT STRUCTURE 

In this section, we investigate the fixed-point struc- 
ture of the RG flow. We first concentrate on the 
local-potential approximation and later include next-to- 
lcading-order terms in the derivative expansion. In fact, 
it turns out that there is a qualitative difference of the 
fixed-point superpotentials between the different orders. 
Similar observations are known from two-dimensional 
bosonic theories (341 [35j and are a particularity of two- 
dimensional systems. Still, the local-potential flow is in- 
teresting in its own right. Its IR flow is also quanti- 
tatively relevant for studies of the phase diagram, see 
Sect. ELD 

Since RG fixed-point studies require a scaling form of 
the flow equation, we switch to dimcnsionlcss quantities 
w and t defined by W k (4>) = kuit{(j>) and t = ln(fc/A). In 
two dimensions, a scalar field is dimensionless such that 
no dimensionful rescaling is required. The flow equation 
(|37| for the dimensionless quantities reads 



dtw t {(j)) + w t {4>) 



1 



(39) 



4tt1 + <((£) 2 

with d t — kd k . The fixed points are characterized by 
dt w '* = 0. In the following, we solve the fixed-point equa- 
tion by various methods. 
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A. Polynomial expansion 



For small values of the field, a polynomial approxima- 
tion for w' t is justified. If w' t ((f>) is an even function at 
the cutoff scale then it remains even at all scales. Its 
expansion reads w' t {4>) — \t{<p 2 — a t) + 64, t0 4 + ^6,t0 6 + 
b s t tfi s + . . . . The dimensionless couplings A t , ba relate to 
the bare couplings A,6^ through A = kX t and 6$ = kb t ^. 
a is dimensionless, therefore we have a = at- 

Expanding both sides of the flow equation in terms 
of (f>, a comparison of coefficients leads to the following 
system of coupled ordinary differential equations: 



d t a\ 

d t Xt 
d t bi.t 
dth, t 



1 

"2^r 



,364,, 
7rA t 



36 4 , t 



6A? 

1 M 

ir ir 

156 6 f 606 4i -A 2 40Af , 
+ ; L 

146 8ii 5606 4 t Xt 16861 At 



(40) 



IT IT 

1266 6 , t A 2 224A[ 



-6 fi 



dtb 2n ,t 



(n + l)(n + 2) 



4tt 



t + /2ra (A*, 64,4, 



Note that only w" enters the right-hand side of the flow 
equation (|39|) : in particular, the lowest coupling constant 
a 2 does not enter the equations for the higher order cou- 
plings. 

At a fixed point, the coupling constants, marked by an 
asterisk, become scale invariant such that the left-hand 
sides in Eq. (|40|) vanish. The corresponding system of 
equations 



4tt 



J 2n+2 



(n+l)(n + 2) 



fin (A*, 64 



,&!„) 



(41) 



can be solved iteratively due to the triangular form of 
the system of flow equations (j40|> . At a fixed point, b* 2n 
is a polynomial of order 2n + 1 in A*. Because of the %i 
symmetry of the system of equations, these polynomials 
are odd. In particular, b\ (A*) = 2A* 3 - 7rA*/3 such that 
the projection of any fixed-point solution on the subspace 
defined by the couplings A* and 6| fall onto the curve 
6|(A*) depicted in Fig. [3] Inserting 6|(A*) into the first 
equation in ((40)) yields (a*) 2 = 1/27T. Later, we shall 
see that a 2 defines an IR unstable direction near a fixed 
point. 

Let us truncate the polynomial expansion of the flow 
equation for wt and keep only terms up to order <fi 2n+1 
in This is equivalent to keeping the lowest n fixed 

point equations which yield b\ m (A*) for m < n. In 
the nth equation, the higher-order coupling 62„+2 occurs 



R-stable fixed-points 



Gaussian fixed-point 





V 



Jo 

bji).2r 3 -itr/3 
0.5 



FIG. 2: Projection of the coefficients of all fixed points for 
different truncations on the plane of the couplings A and 64. 



which we set to zero at the fixed point, b\ nJr2 = 0. This 
leads to the polynomial equation 



/ 2 n(A*) = / 2 „(A*,6:(A*), 



,6^(A*))=0, (42) 



where the solutions b\ m (A*) are to be inserted. With 
Mathematica, we have checked up to order 2n = 22 
that all 2n + 1 roots of the odd polynomial /2 n (A*) of 
order 2n + 1 are real. 

Due to the underlying TLi symmetry, the remain- 
ing fixed points of the truncated system come in pairs 
±(A*, 6|, . . . , 6| n ). Hence, we find n independent non- 
trivial solutions to the fixed-point equations in addition 
to the GauBian fixed point where all couplings vanish. 
As discussed in the next subsection, only one of these 
solutions belongs to an infrared stable fixed point with 
all but one eigenvalues in the stability matrix being pos- 
itive (i.e., all but one critical exponents being negative). 
It turns out that this infrared-stable fixed point corre- 
sponds to the largest root of Eq. (|4"2"1) . as indicated in 
Fig. [5] With increasing order of the polynomial trunca- 
tion the root belonging to the IR-stable fixed point con- 
verges to A cr ;t — 0.9816. Roots belonging to any other 
fixed point are bounded by 

— A cr it < A ' < A cr it- (43) 

In Fig. [3] (left panel), we have plotted the values for the 
two lowest coefficients A* and 6| at the infrared-stable 
fixed point in the polynomial expansion for different or- 
ders of truncations. We observe a rapid convergence with 
increasing order of the polynomial approximation. This 
suggests that the polynomial approximation to the su- 
perpotential in the local potential approximation has ac- 
ceptable convergence properties. 



1 This prescription is not unique. Alternatively, we could set 
&2 ri _i_2' f° r i ns tance, equal to its perturbative one-loop value. In 
any case, the choice &2n+2 = used here is self-consistent in the 
sense that the equations of 6<2n are closed and do not depend 
on further input. 
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FIG. 3: Left panel: The first two coefficients A and i>4 of the IR stable fixed point for different truncations. Right panel: Ratio 
of successive couplings of the IR-stable fixed point yielding an estimate of the radius of convergence. 



Coefficients at IR-fixed point 



2n 


A* 


K 


^6 




Ko 


b 12 


6*4 


2 


0.7236 














4 


0.9019 


0.5227 












6 


0.9535 


0.7354 


0.8372 










8 


0.9711 


0.8148 


1.199 


1.694 








10 


0.9777 


0.8451 


1.345 


2.420 


3.801 






12 


0.9802 


0.8570 


1.402 


2.716 


5.401 


9.030 




14 


0.9812 


0.8617 


1.425 


2.836 


6.054 


12.77 


22.23 



TABLE I: The coefficients of the IR-stable fixed point poten- 
tial for different truncations. 



In Fig. [3] (right panel), we have plotted the inverse ratio 
of successive couplings at the infrared-stable fixed point 
for increasing truncation order. This facilitates an esti- 
mate of the radius of convergence for a series expansion 



of w[ in powers of cj> 



T = lim 

' con -LAiii. 



>2, 



J 2i+2 



An extrapolation leads to the approximate value 



0.35 such that for \<f>\ < r con the polynomial \*((j> 2 — 
2 ) + ^2m^ 2m converges to the IR-stable fixed point 



a 



solution w^. In Table [H the coefficients in the polynomial 
approximations to at the IR-stable fixed point are 
listed. Note that the fixed-point values are generically 
regulator dependent, and thus are not directly related to 
physical quantities. 



1. Stability analysis and critical exponents 

Whereas the values of the fixed-point couplings are reg- 
ulator dependent, the critical exponents are universal and 
give rise to a classification of the fixed points. The criti- 
cal exponents are defined as the negative eigenvalues 8 1 



of the stability matrix at the fixed point 
j _ d(d t bj) 



Br 



b=b* 



(44) 



where we have set bo = a\ , hi = A, and I labels the differ- 
ent critical exponents 9 1 and eigendirections v\. Critical 
exponents with positive real part correspond to RG rel- 
evant directions, whereas exponents with negative real 
part mark irrelevant directions. Inserting the flow of At 
into that of a\ , cf. Eq. (|4"0")) . we obtain 



2tt \ t 



(45) 



We observe that the 00-component of the stability matrix 
at any fixed-point yields, BrP = — 1. This together with 
the fact that the remainder of the first column vanishes, 
Bi>i° = 0, implies that a\ is always an eigendirection of 
B^ with corresponding critical exponent 9° = 1. Note 
that this result is manifestly regulator independent and 
thus universal. We conclude that any fixed point of the 
superpotential in the local-potential approximation has 
at least one RG relevant direction. In analogy with po- 
tential flows near the Wilson-Fisher fixed point of Ising- 
like systems, we introduce the following notion for the 
leading critical exponent corresponding to this relevant 
direction: 



v w 



1 



l. 



(46) 



Even though vy/ plays the same role for the superpoten- 
tial flow as the critical exponent v does for the poten- 
tial flow in Ising-like systems, it should be stressed that 
vw does not correspond to the scaling exponent of the 
correlation length (as v does in Ising-like systems). We 
will later see that vw quantifies certain properties of the 
phase diagram. 

Depending on whether we study the UV or IR flow 
of the system, the fixed points have a different meaning. 
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A* 


Critical exponents 6 1 


±.9816 


-1.54 


-7.43 


-18.3 


—37.3 


-68.9 


— 120 


— 204 


—351 


±.8813 


6.16 


-1.64 


-9.82 


—25.6 


-52.5 


-96.9 


-170 


-300 


±.7131 


21.4 


4.37 


-1.57 


-11.1 


-30.1 


-63.3 


-120 


-223 


±.5152 


28.7 


13.3 


3.33 


-1.39 


-11.6 


-32.8 


-71.7 


-145 


±.3158 


20.0 - 4.55 i 


20.0 + 4.55 i 


8.40 


2.57 


-1.14 


-11.6 


-34.3 


-80.4 


±.1437 


11.2 + 9.02 i 


11.2 - 9.02 i 


8.63 


5.19 


1.95 


-.842 


-11.1 


-35.7 


±.0322 


4.20 + 1.18 i 


4.20 - 1.18 i 


2.86 


2.72 + 6.47 i 


2.72 - 6.47 i 


1.47 


-.540 


-10.5 


±.0003 


1.57 + .125 i 


1.57 - .125 i 


1.43 + .702 i 


1.43 + .702 i 


1.14 


.542 + .982 i 


.542 + .982 i 


-0.221 





1 


1 


1 


1 


1 


1 


1 


1 



TABLE II: Critical exponents 8 1 (negative eigenvalues of the stability matrix) for a polynomial truncation at 2n = 16 for the 
nine different fix points in the local-potential approximation. The first exponent 8° = 1 which is common to all fixed points is 
not shown here. 



Towards the UV, any of the fixed points which we have 
found can be used to define a UV completion of the model 
in the sense of Weinberg's asymptotic safety scenario [361 ]. 
All relevant directions emanating from the fixed point 
span the critical hypersurface. The dimensionality of this 
critical surface, i.e., the number of critical exponents with 
Rc 9 1 > 0, corresponds to the number of physical param- 
eters which have to be fixed in order to unambiguously 
define the flow towards the IR. 2 Once these initial con- 
ditions to the flow are provided, any other quantity or 
correlation function can be predicted within the theory. 
Since 9° = 1, we conclude that any UV completion has 
at least one physical parameter. 

Within the local-potential approximation at order 
2n = 16 of the polynomial expansion, the critical ex- 
ponents 9 1 - 1 are given in Table UT1 for all 17 fixed-point 
potentials. By using a different regulator, we check in 
App. [B] that the regulator dependence of the relevant 
positive critical exponents is rather small (up to 10% or 
much less), which confirms the reliability of the present 
truncation. Classifying the fixed-point potentials by the 
slope A* of the potential w' t as a function of </> 2 at cf) 2 = 0, 
the number of relevant directions increases as the slope 
| A* | decreases. The different fixed-point potentials in the 
local-potential approximation - if they persist to higher 
truncation orders - thus correspond to different UV com- 
pletions of the present system with increasing physical 
parameters. These different UV completions thus define 
different nonperturbatively renormalized Wess-Zumino 
models in two-dimensions. 

As for the flow towards the IR, the fixed points can 
generically be related with critical points in the phase di- 
agram of the system. Since the relevant directions are IR 
repulsive, fine-tuning the relevant direction to the fixed 
point corresponds to tuning the system onto its critical 



point. In this sense, the relevant direction corresponding 
to a 2 with vw = 1/9° = 1 is similar to the temperature 
parameter in Ising-like systems (or a mass parameter in 
0(A^)-type relativistic models). For instance, in the do- 
main of attraction of the maximally IR-stable fixed point 
with only a 2 as relevant direction, the tuning of a 2 dis- 
tinguishes between the supersymmetric and symmetry- 
broken phases of the model. More generally, if a system 
is in the domain of attraction of a fixed-point with N 
relevant directions, the phases of broken and unbroken 
supersymmetry are separated by an V-dimensional hy- 
persurface in the space of couplings. 

There is one important difference to Ising-like systems: 
the coupling a 2 associated with the one common relevant 
direction does not feed back into the flow of the higher- 
order couplings. Therefore, the remaining couplings are 
attracted towards the maximally IR-stable fixed point 
for any regular trajectory irrespective of the flow of a 2 . 3 
We conclude that the maximally IR-stable fixed point 
governs the flow towards the IR of w' t '(<p) in the domain, 
where the polynomial expansion is valid. 

Let us finally mention that the critical exponent vw = 
1/9° = 1 receives corrections at higher orders in the 
derivative expansion, cf. Eq. (|57p . Still the relevance 
of the maximally IR-stable fixed point for the IR flow of 
the potential persists. 

In Tabic IILTI wc collected the eigenvalues (negative crit- 
ical exponents) for the maximally IR-stable fix point for 
different truncations. 



B. Solving of the nonlinear differential equation 

For large values of the scalar field, the polynomial trun- 
cation is not valid anymore. Hence, we consider here 
the full nonlinear ordinary differential equation that de- 



2 For marginal directions with Re 6 1 = 0, the flow in the fixed- 
point regime has to be studied beyond linear order. Depending 
on the sign of the first non- vanishing order, these directions are 
again either marginally relevant or irrelevant. 



Irregular trajectories can run to infinity at a finite value of k. 
These divergences are either physically meaningless or signal the 
breakdown of the truncation. 
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An 


A 


eigenvalues of the stability matrix 


1 


iu.yuiy 










16.35 


1.846 


6 


±0.9535 








42.32 


12.00 


1.716 


8 


±0.9711 






79.83 


30.67 


9.951 


1.635 


10 


±0.9777 




129.1 


58.61 


25.05 


8.794 


1.588 


12 


±0.9802 




190.6 96.49 


47.97 


21.75 


8.101 


1.561 


14 


±0.9812 




264.5 144.8 79.50 


41.51 


19.67 


7.680 


1.546 


16 


±0.9816 


351.2 


204.1 120.3 68.90 


37.25 


18.30 


7.427 


1.539 



TABLE III: Eigenvalues of the stability matrix (negative crit- 
ical exponents) of the maximally IR-stable fix point for dif- 
ferent polynomial truncations in the local-potential approxi- 
mation. The first exponent 6° — 1 is not shown here. For 
the first subleading exponent (last column), the polynomial 
expansion shows a satisfactory convergence. 




-0.4 -0.2 0.2 0.4 0.6 



FIG. 4: Illustration of the discussion of oscillating solutions, 
see text. 



scribes the fixed point potential in the local-potential ap- 
proximation. In the polynomial approximation, the po- 
tential wt and w' t contain the IR-unstable coupling at 
which does not flow into the fixed point. Upon differen- 
tiation of the fixed-point equation for w'^cfi), 



l — w 



tii 



4tt (1 + < 2 ) 



l/2\2 ' 



(47) 



with respect to 0, we arrive at the following fixed-point 
equation for the a* -independent function w"(</>) = u(<j>), 

(1 - u 4 )u" = 2u' 2 (3 - u 2 ) u - (1 + u 2 ) 3 4nu. (48) 

Since a 2 does not appear in Eq. (|48[) . we expect to find 
a fully IR-stable solution to this nonlinear differential 
equation (in addition to further solutions with IR unsta- 
ble directions). 

As before, we consider odd solutions u{4>) which are 
fixed by the initial conditions u(0) = and a finite value 
for u'(0). In fact, we find a continuum of oscillatory 
solutions that are defined for all values of <fi. In addition, 
we identify solutions that hit the singular line u(<f>) = 1 
of the differential equation (|4"8"]) but can be continued 
without cusps. Solutions of the second class exist only in 
a finite 4> range. 



1. Oscillating solutions 

Let us integrate the fixed-point equation (|48]) with ini- 
tial conditions at — <fro < (see Fig.@|: 

u{— 4>o) = uq, u'{—<f>o) — where — 1 < Uo < 0. 

The initial point is an extremum of the solution. The 
differential equation (|48[) implies that this extremum is 
indeed a minimum such that u approaches the <fr axis 
away from — cj>Q. Actually it must intersect the axis at 
some point. Otherwise, it would possess a maximum at 
some (f>i with u(— (j>o) < u { ( t ) i) < or it would mono- 
tonically approach the <p ax is without crossing it. But 



a maximum with — 1 < u((f>i_) < contradicts the dif- 
ferential equation (|4"8|) which implies that u"((f>i) > 0. 
Also, if u approached the cf> axis monotonically from be- 
low then near the axis the differential equation would 
imply u" ~ —Attu > 0. Therefore, the slope of u would 
increase with increasing 4> such that finally u would in- 
tersect the 4> axis. 

We conclude that u must intersect the <fr axis and we 
may assume that this happens at (f> = 0, such that the so- 
lution starts off at the minimum at — (f>Q < and hits the 
axis at the origin. This solution extends to an odd solu- 
tion owing to the symmetry u(<f>) — > — u(— (jj) of the differ- 
ential equation and hence has a maximum at </>o > with 
< u((f)o) = —u(—<fio) < 1. Because of the translational 
invariance — > + c and the symmetry u(<f>) — > u(—<p) of 
the fixed point equation the solution must be symmetric 
relative to i-e. u(4>o — (f>) — u^o +</>)■ This proves that 
every solution with — 1 < u(extrema) < 1 is periodic and 
takes its values between —1 and 1. The lines u(cf>) = ±1 
repel solutions oscillating in the strip — 1 < u < 1 as long 
as the slope u'(0) — 7 is less than 7 crit ~ 1.964. Solutions 
with 7 > 7 cr it hit the singularity at u = 1. 

In a similar fashion, one argues that there exists a sec- 
ond class of solutions of the differential equation having 
just one minimum with u > 1 or having just one maxi- 
mum with u < — 1. Solution in this class are not periodic 
and never hit the singular line u = ±1. Since they cannot 
be continuous and antisymmetric they are discarded. 



2. Comparison with the polynomial expansion 

In the preceding subsection, we have seen in the poly- 
nomial expansions of the truncated system that the slope 
7 = 2A is bounded by 2A cr i t — 1.964. This is just the crit- 
ical value 7 cr it for the existence of oscillating solutions of 
the fixed point equation (|48)l . We conjecture that a poly- 
nomial solution belonging to a fixed point with two or 
more relevant directions of the truncated systems, corre- 
sponding to a non- maximal root of f<z n in J42|) , converges 



11 




FIG. 5: Left panel: Comparison between the numerical solution to the differential equation (ODE) and the polynomial ap- 
proximation (Poly) to 16th order for three different fixed points. The fixed points FP1, FP2, and FP3 have the initial slope 
7 = 0.287, 1.4262 and 1.963. The fixed point FP3 is the maximally IR-stable fixed-point. Right panel: The first derivative of 
the potentials in Fig. [S] 



to the Taylor series of an oscillating solution. In Fig. [5] 
(left panel), we have plotted three full numerical solu- 
tions and the corresponding polynomial approximation 
truncated at </> 16 with the same initial value 7 = 2A. For 
the first half period, we find an excellent agreement be- 
tween polynomial approximation and numerical solution. 



3. Solution with u((f>) = 1 for some field value <j> 

Regular periodic solutions only exist for — Tcrit < 7 < 
7 cr it- Increasing the slope at the origin gradually from 
to 7crit, the value it(</> ma x) at the maximum approaches 
the singular line u = 1 and finally hits the singularity at 
the critical field cr it — 0.3704. At the same time the cur- 
vature at the maximum tends to —00. With increasing 
order the IR stable fixed points for the polynomial trun- 
cations converge to the Taylor expansion of the solution 
with the critical slope 7 cr it- 

In order to study the solutions near the singular line, 
we insert the Taylor expansion u(</> C rit + S(f>) = u(0 C rit) + 
ai5(f) + d2<50 2 /2 + . . . with u((f> cr it) — 1 and compare 
coefficients. One sees that the expansion coefficients are 
finite if a\ = u'(0 cr it) = ±\/87r. If this condition is not 
met, the solution hits the singular line u((f>) = 1 with 
infinite slope, as can be seen by studying the differential 
equation for the inverse function <fi(u). The behavior of 
a solution depends in an essential way on 7: if the initial 
slope is less than 7 cr it then the solution is smooth and 
periodic, if the initial slope is 7 cr it then it hits the singular 
line u(4>) — 1 with slope V87T and if the initial slope 
is bigger than 7 cr it then the solution hits the singular 
line vertically; see Fig. [SI If we viewed the solutions 
hitting the singular line as parametric continuation of 
the periodic solutions as 7 approaches 7 cr it, we would 
reflect the solution at the singular line, similarly to the 
solution FP3 in Fig. [5J left panel. For the maximally IR- 
stable fixed-point solution with initial slope 7„it the slope 



at cr it would then jump from \^8tt to — V&tt, as shown 
in Fig. right panel, where we depicted the function 
u'((f) — w"'(<fi) for three different values of the initial 
slope 7 = 2A. 

However, there is a way to extend the solutions hitting 
the singular line without cusps. To see this more clearly 
we note that v — l/u fulfills almost the identical fixed- 
point equation as u, 

(1 - v 4 )v" = 2v' 2 (3 - v 2 ) v - 4tt(1 + v 2 f/v. (49) 

Upon approaching the singular line u — > 1 and v —> 1, 
the two equations (|4"8|) and (f4"9"|) become identical. This 
implies that near the singular line the reflection at the 
singular line maps solutions into solutions. Thus, all so- 
lutions with 7 > 7 cr it hitting the singular line vertically 
or with slope ^/8tt can be extended without cusps beyond 
the singular line. 

We have studied these solutions for large values u^S> 1. 
It is not difficult to see that there exist no solutions with 
u ~ (f> a for large <j>. (This will become different at next- 
to-leading order in the derivative expansion.) We find 
the asymptotic solution 

Was (0) = e CTrl ( ±2 ^ e ^^ +C2) ) 2 -* (50) 

with ci ~ —20.02 and C2 — —0.423 for the maximally 
IR-stable fixed point extended without cusp beyond the 
singular line. It is finite only for 

- </>crit,2 < < 0crit,2, 0crit,2 = 0.5823. (51) 

An unbounded and cusp-free fixed point solution belongs 
to a field theory with compact target space. 

Let us summarize our findings. As a regular oscillat- 
ing or a non-diffcrcntiablc bouncing solution u((f>) is finite 
for all values of the field, this implies a vanishing dimen- 
sionful W^((f) = ku(4>) as the scale k is lowered to the 
infrared. The solutions which penetrate the singular line 
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andr 2 = (k 2 /p 2 - 1)6(1 -p 2 /k 2 ) 
tities, the flow equations read 



For the rescaled quan- 



dtto' t 



if 



n 



(xKY 



l wi 



An n>'/ 2 



ln(l 



1-21 



to 



l + ml 



r/ In 1 + to'i 



2w? 



112 



l + w' 



2< 



(l + <) : 



(52) 



37/ 

2 



(53) 



FIG. 6: All types of possible solutions to the fixed-point 
differential equation in local-potential approximation: (1) 
7 > 7crit, (2) 7 = 7crit, finite (p range, (3) 7 = 7 cr it, oscil- 
lating solution, (4) 7 < 7crit, oscillating solution, (5) solution 
with just one extremum 



u = 1 without a cusp are unbounded from below and 
above and confine the field to a finite interval. The dif- 
ferent types of solutions are depicted in Fig. [S] 



where we have dropped the arguments in fDt(x) for sim- 
plicity. Since the anomalous dimension is assumed to 
be constant in this approximation, we have projected 
Eq. ([53]) onto \ = 0. Note that the limit tr>"(0) -> 
of the right-hand side of Eq. ([55)) exists, yielding 



AX 2 



77 = 



A 2 + 2tt' 



(54) 



2. Polynomial expansion and superscaling relation 



C. Fixed points at next-to-leading order in the 
derivative expansion 

1. Next-to-leading-order flows 

Two dimensional scalar field theories and their super- 
symmetric extensions exhibit infinitely many scale invari- 
ant fixed-points characterized by the central charge c of a 
conformal field theory (37j . In the local potential approxi- 
mation to scalar field theories, only periodic, sine-Gordon 
type fixed-point solutions are accessible. At next-to- 
leading order in the derivative expansion, one finds ad- 
ditional non-periodic fixed-point solutions [34 , l35l | . Here 
we find analogous results for the two-dimensional Wess- 
Zumino model. 

In this subsection, we consider fixed point solutions 
with scale-dependent but field-independent wave func- 
tion renormalization in the next-to-leading order ap- 
proximation. The corresponding flow equations are de- 
rived in appendix [XJ The full next-to- leading-order ap- 
proximation would include a field-dependent wave func- 
tion renormalization Z/ £ ($). For supersymmetric quan- 
tum mechanics, this order has been computed in [l9l |. 
Here, we confine ourselves to the simpler approximation 
Z fc ($) -> Z fc (0) = z k . 

Let us introduce renormalized fields x, by rescaling <j) 
with the wave function renormalization 4> — > x = Zk4>- 
This implies a dimensionless renormalized superpotential 
WtOt) = Wk{x/Zk)/k. The calculation of the flow is out- 
lined in App. [XI in order to avoid artificial IR singulari- 
ties, different regulator shape functions were used in com- 
parison with the local-potential approximation: ri = 



The polynomial expansion of the next-to-leading-order 
superpotential flow equation in terms of dimensionless 
renormalized couplings is given in App. El For instance, 
the flow of the renormalized parameter a 2 can be written 
as (cf. Eq. ([AT2]) ) 

< W-s( 1 -i)-( 1 -I)'*-s* > " (55) 

which is the next-to- leading-order analogue of Eq. ([45[) . 
The renormalized couplings are related to their unrenor- 
malized analogues by 

Similarly to the local-potential approximation, we ob- 
serve that the 00-component of the stability matrix at 
any fixed point yields Bq° = —(1 — where 77 = rj* has 
to be evaluated at the corresponding fixed point. Since 
the remainder of the first column vanishes, B ( j > >1 = 0, 
the coupling a 2 remains always an eigendirection of 
at any fixed point with a critical exponent 9° = — (1— tt), 
implying a superscaling relation 

where we have again introduced an Ising-like notation 
for the critical exponent of the superpotential associated 
with the a 2 direction. This is a remarkable relation as it 
relates this superpotential exponent with the anomalous 
dimension. Recall that in Ising-like systems the thermo- 
dynamic main exponents (i.e., a, /?, 7 and S) are related 
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2n 


2 


4 


6 


8 


10 


12 


14 


V 

l/v w 


0.3284 
0.8358 


0.4194 
0.7903 


0.4358 
0.7821 


0.4386 
0.7807 


0.4388 
0.7806 


0.4387 
0.78065 


0.4386 
0.7807 



side of the flow equation ([52 



TABLE IV: Numerical verification of the superscaling rela- 
tion (|57[) : anomalous dimension 77 and the critical exponent 
1/uw of a 2 for increasing orders in a polynomial truncation 
evaluated for the maximally IR-stable fixed-point. 



among each other by scaling relations, and can be de- 
duced from the correlation exponents v and r\ by hyper- 
scaling relations. Beyond that there is no general relation 
between v and r\. The superscaling relation (|57[) thus rep- 
resents a special feature of the present supersymmetric 
model. 



We would like to stress that Eq. ([57]) is an exact rela- 
tion to next-to-leading order in the supercovariant deriva- 
tive expansion of the effective action. In particular, the 
inclusion of a field-dependent Zk{<fi) implying 77 — ► rj(4>) 
does not modify the superscaling relation, since the su- 
perscaling relation arises from the expansion in <f> near 
4> = 0. Beyond next-to-leading order, Eq. ([5T]) might, in 
fact, receive corrections, since higher- derivative operators 
can still take influence on the flow of the superpotential 
mediated by higher-order interactions between the scalar 
field and the auxiliary field. Whether or not these in- 
teractions play a role for the superscaling relation at the 
fixed points needs to be clarified by future studies. 

Within the present next-to-leading-order truncation, a 
numerical determination of the critical exponents from 
the full set of polynomially expanded flow equations, of 
course, confirms the superscaling relation to a high accu- 
racy. Numerical values for 77 and vw at the maximally 
IR-stable fixed point for increasing order of polynomial 
truncations are given in Table HVl We also observe a 
rapid convergence of the polynomial expansion, yielding 
our best estimates 77 ~ 0.4386 and l/i>w — 0.7807 for the 
critical exponents at the maximally IR-stable fixed point. 
As the anomalous dimension is comparatively large, we 
expect significant quantitative corrections to arise from 
higher orders in the derivative expansion. 

As discussed below, the superscaling relation has an 
immediate physical consequence for the IR flow of the 
masses in the supersymmetry broken phase. 



Fixed points of the nonlinear superpotential flow at 
next-to-leading order 



In order to go beyond the polynomial expansion, let 
us first study the asymptotic behavior of the right-hand 
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(58) 
(59) 



It turns out to be a self-consistent assumption that these 
asymptotic right-hand sides are subdominant in compar- 
ison with the left-hand side of Eq. ([52]) at a fixed point 
dtto'^ = both for small and large values of \. From this, 
it follows that tt>* is proportional to ^ 2 / r '~ 1 for large to". 
In particular, trj'„ grows faster than any polynomial for 
77 = 0, in complete agreement with our previous results 
in section [VTl 

Now, for a non- vanishing 77 we find a new class of so- 
lutions. For these new solutions, we consider again the 
derivative of Eq. ([52]) . This leads to the following fixed 
point equation for u = tr>" (note that u now contains a 
wave function renormalization in contrast to section fVTl): 



An 



y(3 + u 2 ) 
2u 4 
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= (77- l)u+ — u + — 

Z Z7T 



2 3?7 
h u 2 2u 2 
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77(6 + m 2 ) 
2u 5 
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"75 Ml 



(60) 



For an initial condition in terms of an odd superpotential 
Wa at the UV scale k = A, tt>t is odd and the fixed point 
solution u — to" vanishes at the origin. Thus for weak 
fields we have u(x) ^ 1 an d we may expand the loga- 
rithm in powers of u. The resulting fixed point equation 
for small u is regular for u — > and reads 



^({r 1 -A) + (&-r 1 )u 2 + h -{r 1 -2)u i 



-Xu' - (r?- l)u. 



)u 2 + ^(10- V )u 4 + 



(61) 



Following [34[, we first consider 77 in Eq. ([60[) as a free 
parameter. The initial conditions u(0) = and it'(0) = 
7 = 2A are parameterized by the slope 7 at the origin. 
A solution of Eq. (|60[) with generic slope will run into a 
singularity because the factor multiplying u" eventually 
becomes zero. By fine-tuning the slope it is possible to 
find regular solutions for a given value of r\. In Fig. [7] 
(left panel) we show three regular potentials for rj = 0.1. 
For large values of x, they behave like fo"(x) ~ X 18 - 

These regular solutions define curves of fixed-point so- 
lutions in the 7-77 plane. This is shown in Fig. [7] (right 
panel). For 77 = 2/3 we find a potential that behaves as 
u ~ x m the asymptotic region. We do not find solu- 
tions with larger values for r\. For < 77 < 2/3 it follows 
from simple monotony arguments that the factor multi- 
plying u" in Eq. ([6T)]) has only one node at some value 
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FIG. 7: Left panel: regular potentials for 77 = 0.1. Their asymptotic behavior is to"(x) ~ X 18 i Right panel: Lines of fixed 
points in the 77-7 plane (solid curves) and the anomalous dimension as a function of 7 = 2A obtained from Eq. (|54p (dotted 
curve). We also display the fixed point solutions obtained from a polynomial approximation of Eq. (|52|l and (|53[) for different 
truncations. 



Number of nodes rj vw 

0.4386 1.2809 

1 0.20 1.11 

2 0.12 1.06 



TABLE V: Critical exponents of the first fixed points. 



4>q . The potentials are therefore regular if the right-hand 
side of Eq. (f6U|) vanishes for the same value 4>o , which is 
achieved by a fine-tuning of 7. The outermost curve in 
Fig. [7] (right panel) corresponds to a potential to" with 
no nodes, the next curve to potentials with one node and 
the third curve to potentials with two nodes. We expect 
to find more curves for small 77 and 7 corresponding to 
potentials with more nodes. In Fig. [7] (right panel) we 
also display 77(7) = 47 2 / (7 2 + 87r) obtained from equation 

(53). 

The polynomial approximation to the fixed point solu- 
tion of Eq. (|52p and Eq. ([55)) converges to the maximally 
IR-stable fixed point with 77 = 0.4386 and 7 = 1.759, 
which is just the point of intersection with the line of 
fixed points corresponding to potentials with no nodes. 
The intersection points of the 77(7) line with the other 
lines of fixed points with higher numbers of nodes then 
give estimates for the critical exponents of these other 
fixed points at next-to-leading order in the derivative 
expansion. E.g., we find 77 ~ 0.20 for the fixed point 
with one node and 77 ~ 0.12 for the fixed point with two 
nodes. The corresponding critical exponent i>\y then fol- 
lows directly from the superscaling relation (|57[) . They 
are listed in Tab. |V] The point characterized by rj = 
and 7 = 3.529 in Fig. [7] (right panel) (see appendix IB"]) 
belongs to a solution of the type discussed in section fVTl 
where the maximally IR-stable solution can be extended 
without cusps beyond the critical value of the potential 



D. Synthesis: derivative expansion results at 
leading order and next-to-leading order 

At a first glance, the fixed-point potentials obtained 
at the various orders in the derivative expansion seem 
even qualitatively different. At leading-order, we find 
oscillating solutions, solutions with cusps, and solutions 
with a compact target space. By contrast, the next-to- 
leading-order solutions admit superpotentials that can 
be extended to infinite field amplitude with a standard 
powerlaw asymptotics W'{(j) — * 00 ) ~ </> 2 / ?J ~ 1 — > 00. 

Of course, there are also many similarities, as the 
next-to-leading-order fixed-point potentials can be clas- 
sified by their number of nodes, i.e., they typically ex- 
hibit an oscillating behavior for small fields. In addi- 
tion, we expect the occurrence of superpotentials with 
singular structures at finite field values which have not 
been searched for as systematically as in the leading- 
order case. 

The key to a unified understanding of both orders is 
provided by the anomalous dimension 77, as the new large- 
field asymptotics at next-to-leading order is induced by 
a nonzero value for 77. It should be kept in mind that we 
deduce the nonzero value for 77 from a small-field expan- 
sion of the full flow of the wave function renormalization 
Zk{4>). Beyond this expansion, the anomalous dimen- 
sion will acquire a field dependence 77 — > ij(4>). From 
the general form of the flow equation, we expect that 
the large-field limit is characterized by rj((f> — > 00) — ► 0. 
We therefore conjecture that the true large-field asymp- 
totics of the fixed-point superpotentials lies in-between 
the leading- and next-to-leading-order results. More pre- 
cisely, we expect that a standard asymptotic behavior 
W'(4> — > 00) — > 00 persists, but the powerlaw behavior 
(jp/v—t m ay be replaced by a stronger divergence. 
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In any case, the rapid convergence of the polynomial 
expansion in both orders of the derivative expansion, as 
well as quantitative agreement between observables de- 
rived from the polynomial expansion and from the full so- 
lution support the reliability of the overall picture arising 
from the derivative expansion. 



VII. THE GAUSSIAN WESS-ZUMINO MODEL 

In principle, each of the fixed points defines a different 
UV completion of the Wcss-Zumino model and therefore 
a different physical system with a different number of 
physical parameters. In the following, we concentrate 
on the Wess-Zumino model defined at the Gaufiian fixed 
point corresponding to an asymptotically free theory. At 
least seemingly, this is a natural choice, as it has also 
often been used in lattice computations. However, the 
Gaufiian fixed point actually has infinitely many relevant 
directions, as is already revealed by perturbative power- 
counting. As a consequence, there are strictly speaking 
infinitely many physical parameters. 

In practice, one usually starts with a classical superpo- 
tential including quadratic perturbations of the Gaufiian 
fixed point, W' K = Aa(</> 2 — a\), at the UV cutoff k = A, 
implying that infinitely many couplings have been set 
to zero at that scale. Since the RG trajectories are not 
regulator independent, it will still be difficult to com- 
pare our results with those of, say, lattice computations, 
as the same physical system with lattice regularization 
might have a very different action at the lattice cutoff 
A = n/a. In fact, we find a substantial quantitative 
regulator dependence for non-universal quantities within 
the functional RG calculations; see App. 151 which might 
imply that a meaningful comparison with other methods 
should only be made on a qualitative level. 

On the other hand, one may interpret the choice of 
the regulator as belonging to the definition of the theory 
itself: the regulator together with the initial condition 
in the form of a quadratic perturbation specifies the RG 
trajectory uniquely also at a finite scale A. In the general 
case, this viewpoint has the disadvantage that a change 
of the cutoff scale A on the line of constant physics generi- 
cally involves an adjustment of the couplings of infinitely 
many operators. In order to find out how much these 
operators actually affect the flow of the superpotential, 
we have varied the cutoff scale A. For a given A, we ad- 
just the couplings in W' h = \\{(j) 2 — a\) such that we 
obtain fixed reference couplings a\ and Aa at a refer- 
ence scale Ao, ignoring higher-order couplings. This way 
of looking at the cutoff-dependence is very much moti- 
vated by similar procedures in recent lattice simulations 
of the two-dimensional Wess-Zumino model in [J, |38| . 
For large enough A, the solutions of the flow equation 
show that the dependence of, e.g., the ground state en- 
ergy at k — on the actual cutoff scale is small. This 
observation helps making the viewpoint of including the 
regulator in the definition of the theory practicably ap- 



plicable. Still, it has to be emphasized that we observe 
a significant regulator dependence of the nonuniversal 
quantities, see App. [BJ 



A. Numerical solution of the flow equation in 
local-potential approximation 

At the cutoff scale A, we start with Wj^ = Aa(</> 2 — a\), 
where Aa0 2 is a relevant perturbation of the Gaufiian 
fixed point. For large values of the field, the right-hand 
side of Eq. p8|) vanishes such that d^W^ ~ and the ef- 
fective potential remains close to the classical potential. 
It follows that the dimensionless potential wt = Wk/k 
diverges for large <f> in the infrared. On the other hand, 
we expect that the dimensionful W"{4>) — ku t ((j)) con- 
verges to zero for small fields, since Ut can be attracted 
by the maximally IR-stable fixed point and converges to 
a bounded function u. 

Because of the singularity in the fixed point equation 
for u, the solution of the partial differential equation ([38)) 
poses a numerical challenge. To meet this challenge, 
we exploit the fact that the polynomial solutions with 
scale-dependent coefficients (pR))) yield excellent approx- 
imations to the solution of the full partial differential 
equation for small fields <f> < 4> C Ht- For the fixed-point 
solution, this has been demonstrated earlier, see Fig. 
left panel. Thus, we use a polynomial approximation 
for \<j>\ < 0crit- For large fields \4>\ > C rit, the polyno- 
mial approximation fails and we solve the partial differ- 
ential equation numerically. As boundary conditions for 
the numerical solutions, we impose Wj£(</>) = W^{<f) for 
\<j)\ — > oo and Wk{4>mt) equal to the polynomial approx- 
imation at cr ;t and scale k. Figure [5] (left panel) shows 
the flow of such a potential. 



B. The phase diagram 

Depending on the parameters in Wi — \h{4> 2 — 5a.) we 
may end up with a broken or unbroken supersymmetry 
in the infrared. Note that Afc is dimensionful, whereas afc 
is dimensionless. The rescaled dimensionless couplings 
have an index t = lnfe/A, e.g. A& = k\ t and afc = a t . 
In the following, we determine the parameter region for 
which supersymmetry is dynamically broken. A criterion 
for supersymmetry breaking is provided by a nonvanish- 
ing ground state energy, given by the minimal value of 
V{4>) = \W' 2 (4>), where W = W^ Q . The ground state 
energy is nonzero if and only if W'((f>) > for all <p. Since 
W'(<j)) is minimal at the origin we may use the polynomial 
approximation for which W^(Q) = — Afc • a 2 = —(kXt) ■ a 2 . 
The dimensionless coupling A* flows to the IR fixed-point 
value A*, whereas the coupling a 2 diverges for k — ► or 
equivalently for t — ► — oo, such that the dimensionful 
quantity k A t - a 2 converges to a finite value. This is a 
direct consequence of the fact that the divergence of a 2 
at the fixed point is governed by the critical exponent 
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FIG. 8: Left panel: Flow of a potential with dynamical supersymmetry breaking, W^(4>) = Xa(4> 2 — «a)> A = 1000, Aa = 100, 
a A = 0.3; Right panel: Phase diagram in the space of the dimensionless couplings specified at the cutoff scale A for different 
truncations. 



6 a = 1. The supersymmetric phase is characterized by 
af — > +oo and a double well potential W k 2 , whereas the 
broken phase is characterized by a 2 —* — oo and a single 
well potential. 

This gives rise to a strong analogy between dynami- 
cal supersymmetry breaking in this and many other su- 
persymmetric systems and quantum critical phenomena 
in strongly correlated fermion systems [l|j]: both phe- 
nomena are governed by a control parameter of the (su- 
per-)potential which plays the role of a bosonic mass 
term; see, e.g., [39( for an RG treatment of a semi-metal- 
superfluid quantum phase transition in a fermionic sys- 
tem. This control parameter, which is usually called 8 in 
quantum critical phenomena, is associated to the combi- 
nation S = Aa ■ « A /A in our case. The critical value <5 cr 
of the control parameter marks a quantum critical point 
of a quantum phase transition. 

In Fig. [8] (right panel), we depict the phase diagram 
in the space of dimensionful couplings Aa and Aa • a\ in 
units of the UV cutoff A for increasing truncation order. 

In the (j) 2 truncation, the system of differential equa- 
tions reads (see equation (|40|) ) 



dtaf = — - 



1 

2^ 



6A? 



dtXt = 



■At + ^t, (62) 

7T 



which can be solved analytically. The phase transition 
curve in the (Aa<x a , A) plane is given by the initial values 
for which a 2 ^_ oa changes sign. This condition yields 



Aa&a arcsin(a) 
A ~ x/24^a ' 



7rA 

6Aa' 



(63) 



In order to find the approximate phase transition curve 
for the higher order polynomial truncations, we have in- 
tegrated the corresponding systems of flow equations nu- 
merically. From the (f> 2 to the </> 4 truncation, the transi- 
tion curve moves considerably upwards because the cou- 
pling &4 t enters the flow equation for a 2 in Eq. (|4"D")) . 



The higher couplings &6,t) &8,tj • • ■ enter the differential 
equation for at only indirectly, and, as a result, the ap- 
proximate phase transition curves converge rapidly with 
increasing order of the truncation. We find that there 
exists a critical value for the cutoff parameter AAa A | cr it 
characterizing the phase transition in the strong coupling 
limit Aa — > oo. To lowest order, Eq. ([63]) leads to a crit- 
ical value of AAa A | cr it/A = ~ 0.181. Our best 
estimate for this critical value derived from a numerical 
higher-order solution is AAa A | C rit/A ~ 0.263. We con- 
clude that supersymmetry can never be broken dynami- 
cally above this critical value. This agrees qualitatively 
with earlier results in the literature [23, l2a| . 



A more quantitative comparison to lattice simulation 
is inflicted by the strong regulator dependence of nonuni- 
versal quantities, such as the bare critical coupling val- 
ues discussed above, see also App. |Bj For instance in 
[241 [25|], the phase diagram and the ground state en- 
ergy of the present model were investigated anew. For a 
fixed Aa/A = 0.5, a phase transition from a state with 
broken to a state with unbroken supersymmetry was ob- 
served at the critical value (Aa 2 )| cr it/A = 0.48. With 
two different methods they obtained 0.40 and 0.52 for 
this critical value but they state 0.48 to be the most 
solid value. In these works, the thermodynamic limit has 
been performed without an accompanying continuum ex- 
trapolation. On the other hand, the quoted numbers for 
(Xa 2 )\ were in reasonable good agreement with earlier 
results obtained with the help of the worldline path in- 
tegral method [28]. As demonstrated in App. [Bl a direct 
com par ison of bare quantities between the present work 
and |24l . [2o| is anyway not meaningful due to the scheme 
dependence. If more lattice points in the phase diagram 
were available, dimensionless coupling ratios could be 
compared which are likely to be less affected by scheme 
dependencies. 
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C. Mass scaling 

The scale-dependent bosonic mass can be read off from 
the bosonic potential in the on-shell formulation, which 
is given by Vk{4>) = \{W k {<f))) 2 in our truncation. In 
terms of renormalized fields x = Zk<t>, the dimensionful 
renormalized bosonic potential and bosonic mass thus is 

Vk(x) = l^(W!,(x/Z k )) 2 , ml = Vi'(Xnna), (64) 

where Xmin denotes the minimum of the effective poten- 
tial Vfc(x)' The true mass can then be read off in the 
limit to = lim/c^o mk . 

In the broken phase, the minimum of both 14 and W' k 
is at x = 0, such that the bosonic mass yields 

ml 4^(«"(0) = 2fc 2 A t 2 |a 2 |, (65) 

where we have used the fact that the renormalized pa- 
rameter a 2 is negative. Assuming that the system is 
dominated by the maximally IR-stable fixed point with 
At — ► A* and a 2 ~ k~ x l Vw , the renormalized mass scales 
as 

to 2 - k 1+ i (66) 

where we have employed the superscaling relation (|57jl. 
For 77 > — 2, the renormalized bosonic mass scales to zero 
upon attraction of the maximally IR-stable fixed point. 
Indeed, for the Gaufiian Wess-Zumino model considered 
here, we observe that the flow in the broken phase is al- 
ways attracted by the maximally IR-stable fixed point. 
Together with the fact that the broken phase also goes 
along with a massless goldstino, we conclude that the 
broken phase remains massless in both degrees of free- 
dom. 

In a certain sense, the underlying limit k — > repre- 
sents an extreme point of view. Any experiment as well 
as any lattice simulation will involve an IR cutoff scale fc m 
characterizing the measurement, e.g., the scale of a mo- 
mentum transfer, the detector size or the lattice volume. 
Any measurement therefore is not sensitive to k — > but 
to k — > fc m > 0. We conclude that any measurement of 
the bosonic mass in the broken phase will give a nonzero 
answer proportional to the measurement scale, whereas 
the goldstino will be truly massless. 

In the broken phase, the superscaling relation also has 
a special consequence for the unrenormalizcd potential. 
We observe that 

W'(0) = -X k at = -kZ k \ t a 2 t ~ fc 1 ^ -> const., 

(67) 

where we have used Eq. ([57]) and Z k ~ k n ' 2 . There- 
fore, the flow of the superpotential freezes out near the 
origin if the system is in the domain of attraction of the 
maximally IR-stable fixed point. 
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FIG. 9: Renormalized mass at different scales k as a function 
of the initial condition XacLa a t an initial coupling of A a = 0.1. 



In the supersymmetric phase, the bosonic as well as 
the fermionic mass is given by 

ml = i (WHiXmin/Zk)) 2 , (68) 

where we made use of the fact that W k vanishes at the 
minimum. For typical flows, W' k ' remains positive at the 
minimum also in the k — * limit. This leads to a generic 
decoupling of the massive modes, once k drops below 
the mass threshold. For an accurate inclusion of this 
mass decoupling, also the anomalous dimension needs to 
be evaluated at the minimum Xmim implying rj — > for 
k 2 <m£, cf. Eq. ((531). 

For simplicity, we confine ourselves to the local- 
potential approximation and set Z k — ► 1, f] — ► for a 
computation of the masses across the phase transition. 
Our results are displayed in Fig. [H] as a function of the 
initial-condition parameter Aa^a °f the relevant direc- 
tion. The calculation has been performed at an initial 
value for the coupling of Aa/A = 0.1. We observe a 
critical value of AaS 2 | C rit/A ~ 0.045, above which super- 
symmetry remains unbroken and the theory is massive. 
Both masses vanish at the quantum phase transition. Be- 
low the critical value, supersymmetry is broken and the 
goldstino is massless. The bosonic mass also approaches 
zero for k — ► 0, but remains finite for any finite value of 
k, potentially representing a measurement scale. 

VIII. CONCLUSIONS 

In this work, we have constructed a manifestly super- 
symmetric flow equation for the application of the func- 
tional RG to the two dimensional Af = 1 Wess-Zumino 
model. The regularization turns out to be similar to the 
one in supersymmetric quantum mechanics. This is not 
surprising since supersymmetric quantum mechanics can 
be derived from the Wess-Zumino model by a dimen- 
sional reduction. 
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For approximate solutions to the flow equation, we 
have employed an expansion of the effective action in 
terms of field operators containing increasing powers 
of supercovariant derivatives. This provides a system- 
atic approximation scheme that preserves supersymmc- 
try. We obtain a flow equation for the superpotential and, 
at next-to-leading order, for the wave function renormal- 
ization. 

The Wess-Zumino model has a highly nontrivial fixed- 
point structure in terms of scaling solutions for the su- 
perpotential. At leading- as well as next-to-leading order 
in the derivative expansion, we find fixed-point superpo- 
tcntials which can be classified by an increasing number 
of relevant directions. At leading order, the classification 
of fixed-point solutions can largely be done on analyti- 
cal grounds. Here, we find oscillating sine-Gordon-type 
or even bouncing solutions on the one hand, and solu- 
tions confining the field to a compact target space on 
the other hand. For the regular oscillating solutions, the 
large-field asymptotics is turned into a standard form 
W'(4> — > oo) — > oo at next-to-leading order, where the 
fixed-point superpotentials can be classified also accord- 
ing to their number of nodes. This is reminiscent to 
fixed points of the effective potential in two-dimensional 
bosonic theories which can be related to conformal field 
theories [H], HH| • Exploring the connection between the 
present supersymmetric models at their fixed points and 
conformal field theories remains an interesting question 
for future work. 

Each fixed point defines its own universality class, the 
physics of which is determined by the RG relevant direc- 
tions of the fixed points. In other words, each fixed point 
defines a different Wess-Zumino model with the number 
of physical parameters given by the number of relevant 
directions. One extreme is provided by the Gaufiian fixed 
point which has infinitely many relevant directions in 
agreement with perturbative power-counting arguments. 
The other extreme is given by the maximally IR-stable 
fixed point which has only one relevant direction. We 
have demonstrated that this relevant direction is exactly 
given by a 2 and is shared by all other fixed points as well. 

At leading- and ncxt-to-lcading order, the critical ex- 
ponent associated with this relevant af direction can be 
computed exactly. At next-to-leading order, this yields 
an intriguing relation between this critical exponent of 
the superpotential vw = 1/6° and the anomalous di- 
mension rj of the field. As such a relation between criti- 
cal exponents associated with correlation functions is not 
known from Ising-type systems (where thermodynamical 
exponents are related by scaling relations among each 
other and by hyperscaling relations to the correlation 
exponents), this superseding relation appears to be a 
unique property of supersymmetric theories. 

The superscaling relation has a direct consequence for 
the renormalized superpotential W at the origin in field 
space. It dictates a freeze out of the derivative of the 
superpotential at the origin in the deep IR, if the system 
is governed by one of the fixed points. 



As an example for a model with dynamical supersym- 
metry breaking, we considered the Wess-Zumino model 
defined by a quadratic perturbation of the Gaufiian fixed 
point. In addition to the initial coupling value A|a, the 
model has a control parameter Ao 2 |a the value of which 
decides about the realization of supersymmetry by the 
ground state of the theory. At a given coupling, super- 
symmetry can only be broken for Ao 2 |a below a certain 
critical value which is in accord with a general argument 
by Witten. We have computed the critical line of quan- 
tum phase transitions in the coupling-control-parameter 
plane for a wide range from weak to strong coupling. 
Most importantly, the control parameter stays finite even 
at arbitrarily large coupling. 

We have also computed the masses of the lowest 
fcrmionic and bosonic excitations across the quantum 
phase transition. In the supersymmetric phase, both 
masses are equal and nonzero, but drop to zero at the 
phase transition. In the broken phase, the fermion has 
a massless goldstino mode. We observe that the boson 
also becomes massless in the broken phase in the deep 
IR, k — > 0, but stays finite for any finite k. Associating 
k with a typical measurement momentum scale (say, in- 
verse length of a detector), our results predict that the 
bosonic mass in the broken phase is proportional to the 
momentum scale set by the detector. 

The critical properties of the quantum phase transition 
remain an interesting open problem. In this work, we 
have considered the ground state energy or the bosonic 
and fermionic masses of the theory as order parameters 
for the symmetry. In our numerical results, we have 
found no hint for typical scaling behavior of these quan- 
tities near the phase transition so far. Actually, a true 
field-valued order parameter is given in terms of the ex- 
pectation value of the auxiliary field F. It is therefore 
natural to expect that order-parameter fluctuations of 
the F field play an important role near criticality and 
eventually establish a scaling behavior. Technically, this 
requires the inclusion of potential terms for the F field. 
As these appear at higher orders in the supercovariant 
derivative expansion, a quantitative description of the 
critical regime remains a technical challenge. 

From the perspective of the functional-RG tool box, 
we have solved the flow equation for the superpotential 
both in a polynomial expansion as well as with a full nu- 
merical solution of the corresponding partial differential 
equation. Whereas the polynomial approximation for the 
potential is only reliable in the vicinity of its expansion 
point with a finite radius of convergence, it often suffices 
to extract reliable quantitative information about phys- 
ical observables such as critical exponents or the phase 
diagram. 

In the context of the phase diagram of the Gaufiian 
Wess-Zumino model, we have expressed our concern that 
a quantitative comparison with other methods such as 
lattice simulations can be plagued by the fact that the 
model has infinitely many relevant directions. As a di- 
rect consequence, the results of the model defined by a 
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certain (say quadratic) perturbation of the fixed point at 
a fixed UV scale A are regulator dependent, as this defin- 
ing initial condition is not universal. A much better com- 
parison could arise from defining a Wess-Zumino model, 
e.g., in the vicinity of the next-to-maximally IR-stable 
fixed point, where there are only two relevant directions 
and thus two tunable physical parameters. Fixing these 
parameters in terms of two observables, all other quanti- 
ties are a universal scheme-independent prediction of the 
theory. 
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APPENDIX A: FLOW EQUATION WITH WAVE 
FUNCTION RENORMALIZATION 



At next-to-leading order in the derivative expansion, 
a field-independent wave function renormalization is in- 
cluded in the truncation via 

T k [<l>,F,^,i>] = J ' d 2 x{Z 2 (\d^<\>+\i>$^-\F 2 ) 
+ \W{!{^ 1 ^-W' k {4>)F]. (Al) 
The cutoff action reads 



AS, - I {cj>,F)ZtRl dx^Z'iR^, (A2) 



1 



with Rf and R% given in Eq. (f28|) . The flow equation 
for the superpotential is obtained by a projection onto 
the terms linear in the auxiliary field and by integration 
with respect to <j). We obtain a similar flow equation as 
in Eq. dHSJ), 



1 fd 2 p {l + r 2 )Zld k { ri ZD 



1 fd 2 P (W^^ + nZpd^Z 2 ) 

2 / 4tt 2 A 



(A3) 



Including the wave function renormalization, the expres- 
sion in the denominator reads 



A = z£p 2 (l + r 2 ) 2 + (W{ + r 1 Z 2 ) 2 



(A4) 



For the flow of the wave function renormalization, we 
project the flow equation onto the terms quadratic in the 
auxiliary field. As we consider only a field-independent 



wave function renormalization, we can also project onto 
(b = 



d k Z* = -W' h \$fZlj 0(l + r 2 )x 
^Zl{Wl!{4>)+r 1 Zl){l + r 2 



A 3 



dkinzi) 



(A5) 



zyji + r^-jw^^ + nziy 

A 3 



dk(r 2 Z 2 k ) 



t,=o 



As regulator shape functions, we choose r\ = and r 2 — 
(k 2 /p 2 — 1)9(1— p 2 1 k 2 ) for which the momentum integrals 
in the flow equations can be calculated analytically. We 
obtain 



d k W k 



Wfd k (k 2 Z 2 k ) + k*Z£d k Z 2 



8nW, 



1,3 



In 1 



d k Z 2 k 



4tt 



k 2 d k Z 2 

( z k wf 

W" 2 



Wfkd k Z 2 



-kd k Z 2 k \n 1 



wt 



k*Z* 




™^Z 2 



(Wf 



k 2 Zlf 



J th=0 



For the renormalized fields, \ — Z k (f>, the superpo- 
tential scales as W k (x) = W k (4>) and W k ' = W' k /Z k , 
W k " = Wg/Zl, ... In terms of the anomalous dimension 
rj = —d t Z 2 /Z kl the preceding flow equations read 

kd k W k (x) =lxK - ^777 (A8) 



(77 - 2)k 2 Wf + r)k A 



8tt#; 
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In terms of the dimensionless superpotential lf(x) 
W(x)/k, this reads 



d t tt>k{x) 



K) k 



V 



(i-2K 2 



In 1 




(A10) 
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(All) 



J <*=o 



which agrees with Eqs. (|S"2")) and (|53^) . 
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1. Polynomial approximation 

Next, we perform a polynomial expansion of the super- 
potential flow Eq. (IA10|) including wave function renor- 
malization. We use the conventions to' t (x) = X t (x 2 — 



) + J2n=i b2n,tX 2n i leading to the following system of 



coupled equations: 



. (4->ft)A? 
2tt 



(A12) 



dtdt =^-(4 - m) ~ weft 



8tt 
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Again, the first-order coupling at does not influence the 
higher-order flows. Setting the left-hand side to zero, 
we find a nonlinear system of algebraic equations for the 
fixed-point couplings. The solutions determine the coeffi- 
cients of the fixed-point superpotential in the polynomial 
expansion. 

For a truncation at 2nth order, we again find 2n + 1 
real fixed points. As proved in the main text, each fixed 
point has at least one relevant direction which is provided 
by the a 2 direction. The maximally IR-stable fixed point 
has only this relevant (IR-unstable) direction. The other 
fixed points can be classified according to their increasing 
number of further relevant directions. The Gaufiian fixed 
point remains fully UV attractive also at this order of the 
derivative expansion. 



APPENDIX B: REGULATOR DEPENDENCE 

In this appendix, we repeat several studies of the 
local-potential approximation as done in the main text 
but now for a different regulator shape function T2 — 
(k 2 /p 2 — l)9(p 2 /k 2 — 1). This is actually equivalent to the 
regulator used for at next-to-leading order in the main 
text. On the one hand, this regulator comparison gives a 
rough estimate of how nonuniversal quantities may vary. 



coefficients at IR-fixed point 
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A* 




1,' 




b 10 


b 12 


b* 4 


2 


1.023 
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1.405 


1.301 
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1.540 


2.040 


3.097 
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1.593 


2.374 


4.868 


8.651 








10 


1.615 


2.523 


5.725 


13.38 


26.00 






12 


1.625 


2.590 


6.124 


15.70 


39.55 


81.19 




11 


1.629 


2.620 


6.308 


16.79 


46.11 


121.8 


259.4 



TABLE VI: Coefficients of the first few couplings in a polyno- 
mial expansion of the fixed-point superpotential at the max- 
imally IR stable fixed-point for different truncations. The 
expected nonuniversal deviations from Tab. [I] are due to the 
use of a different regulator. 



On the other hand, it supports all structural results of 
the main text, which should in any case be universal. 

The flow equation for the superpotential ()36|) now 
takes the form 



In 



(Bl) 



In terms of the dimensionless quantity kw t = Wk, this 
reads 

ln(l + <(0) 2 ) 

= - • (B2) 

This equation is regular at w" — since ln(l + x) = 
x — ^- + ^- =F • • • for small x. The second derivative of 
this equation at the fixed point dtw* = is given by 
(« = <) 



Attu 6 (u 2 + 1) 



(u 2 + l)ln {u 2 + 1) - 2u 2 

2u (u 2 - 1) u' 2 



(B3) 



2u'- 



{u 2 + 1) (2u 2 - In (m 2 + 1) {u 2 + 1)) 



which agrees with Eq. (foT))) in the limit r\ = 0. It has the 
same singularity structure as the corresponding Eq. ([48]) 
for the simpler regulator used in the main text. The 
fourth derivative becomes singular if 



1) ln(u 2 



1) 



1.9803, (B4) 



(to be compared with the singular point u = 1 for the 
simpler regulator). Overall, we find the same types of 
solutions as discussed in section I VII In particular, we 
again find oscillatory solutions and a maximally IR-stable 
fixed point with one relevant direction. 

A polynomial expansion of Eq. (|B3[) with u = 2X 2 (f> + 



S^=2 ^ n ' ^2n,*0 2n_1 yields a system of algebraic equa- 
tions that determine the expansion coefficients of the 
fixed-point potential. This system is given by Eq. (|A13[) 
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FIG. 10: Left panel: fixed-point superpotentials calculated from Eq. (|B3|) with A = 1.6295 for the maximally IR-stable fixed- 
point and A = 1.0411 for an IR-unstable fixed-point. Right panel: Comparison between the critical bare values for (\a 2 )t=o 
characterizing the phase transition calculated with the regulators T2 = (k/p— l)6(p 2 /k 2 — 1) and T2 = (k 2 /p 2 — l)6(p 2 /k 2 — 1). 
The difference is attributed to a strong scheme dependence of this nonuniversal quantity. 
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TABLE VII: Critical exponents 8 1 (negative eigenvalues of the stability matrix) for a polynomial truncation at 2n = 16 for the 
nine different fix points in the local-potential approximation. The first exponent 0° — 1 which is common to all fixed points 
is not shown here. This table should be compared with Tab. [TT] which has been computed with a different regulator. For all 
positive critical exponents, we find a remarkable degree of universality, as these exponents for the different regulators differ if 
at all by at most 10%. The subleading negative critical exponents show larger variations. 



without the equation for a 2 , with r\ — and the deriva- 
tives with respect to t set equal to zero. 

The results for the fixed-point couplings for the maxi- 
mally IR-stable fixed-point are displayed in Tab. I VII for 
different truncations. A comparison with Tab. U reveals 
the regulator-dependent variations of these nonuniversal 
quantities. 

We find that the coupling A* converges quickly. Again, 
the polynomial expansion provides a good approximation 
to the solution of the partial differential equation for the 
first half period. This is displayed in Fig. left panel. 

By contrast, the universal critical exponents B 1 at the 
fixed points are much less regulator dependent. This is 
demonstrated in Tab. I VIII which should be read side by 
side with Tab.HTlfor a different regulator. For all positive 



critical exponents, we find a remarkable degree of uni- 
versality, as these exponents for the different regulators 
differ - if at all - by at most 10%. The subleading neg- 
ative critical exponents show larger variations and thus 
require higher orders in the derivative expansion for a 
quantitatively reliable prediction. 

We also calculate the phase diagram with this regula- 
tor for a truncation at 2n = 10 and compare both reg- 
ulators r 2 — (k/ \p\ — l)9(p 2 /k 2 — 1) and r 2 = (k 2 /p 2 — 
\)9(p 2 /k 2 — 1) in Fig.[Tni right panel. The nonuniversal 
values for Xa at the phase transition differ roughly by 
a factor of two. This clearly demonstrates that a naive 
comparison of bare couplings is substantially inflicted by 
the regularization scheme. 
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